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ABSTRACT 

The explanation of the accelerated expansion of the Universe poses one of the most 
fundamental questions in physics and cosmology today. If the acceleration is driven 
by some form of dark energy, and in the absence of a well-based theory to interpret 
the observations, one can try to constrain the parameters describing the kinematical 
state of the universe using a cosmographic approach, which is fundamental in that it 
requires only a minimal set of assumptions, namely to specify the metric, and it does 
not rely on the dynamical equations for gravity. Our high-redshift analysis allows us 
to put constraints on the cosmographic expansion up to the fifth order. It is based 
on the Union2 Type la Supernovae (SNIa) data set, the Hubble diagram constructed 
from some Gamma Ray Bursts luminosity distance indicators, and gaussian priors on 
the distance from the Baryon Acoustic Oscillations (BAO), and the Hubble constant 
h (these priors have been included in order to help break the degeneracies among 
model parameters). To perform our statistical analysis and to explore the probability 
distributions of the cosmographic parameters we use the Markov Chain Monte Carlo 
Method (MCMC). We finally investigate implications of our results for the dark en- 
ergy, in particular, we focus on the parametrization of the dark energy equation of state 
(EOS). Actually, a possibility to investigate the nature of dark energy lies in measuring 
the dark energy equation of state, w, and its time (or redshift) dependence at high ac- 
curacy. However, since w{z) is not directly accessible to measurement, reconstruction 
methods are needed to extract it reliably from observations. Here we investigate differ- 
ent models of dark energy, described through several parametrizations of the equation 
of state, by comparing the cosmographic and the EOS series. The main results are: 
a) even if relying on a mathematical approximate assumption such as the scale factor 
series expansion in terms of time, cosmography can be extremely useful in assessing 
dynamical properties of the Universe; b) the deceleration parameter clearly confirms 
the present acceleration phase; c) the MCMC method provides stronger constraints 
for parameter estimation, in particular for higher order cosmographic parameters (the 
jerk and the snap), with respect to those presented in the literature; d) both the esti- 
mation of the jerk and the DE parameters, reflect the possibility of a deviation from 
the ACDM cosmological model; e) there are indications that the dark energy equa- 
tion of state is evolving for all the parametrizations that we considered; /) the q(z) 
reconstruction provided by our cosmographic analysis allows a transient acceleration. 

Key words: Gamma Rays : bursts - Cosmology : distance scale - Cosmology : cos- 
mological parameters 

1 INTRODUCTION 

At the end of the '90s observations of high redshift supernovae of type la (SNIa) revealed that the universe is now expanding 
at an accelerated rate. This surprising result has been independently confirmed by observations of small scale temperature 
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anisotropies of the cosmic microwave background radiation (CMB) (Riess et al. 2007; Astier ct al. 2006; Kowalski ct al. 2008 
|Spergel et al. 200~7j ). It is usually assumed that the observed accelerated expansion is caused by a so called dark energy, 
with unusual properties. The pressure of dark energy pd e is negative and it is related to the positive energy density of dark 
energy ede by pde = w ede where the proportionality coefficient w < 0. According to the present day estimates, about 75% of 
matter-energy in the universe is in the form of dark energy, so that now the dark energy is the dominating component in the 
universe. The nature of dark energy is not known. Proposed so far models of dark energy can be divided, at least, into three 
groups: a) a non zero cosmological constant, in this case w — — 1, b) a potential energy of some not yet discovered scalar 
field, or c) effects connected with non homogeneous distribution of matter and averaging procedures. In most scenarios of 
dark energy, its properties are mainly characterized by the equation of state parameter (EOS), w. Extracting the information 
on EOS of dark energy from observational data is then at the same time a fundamental problem and a challenging task. For 
probing the dynamical evolution of dark energy, under such circumstance, one can parameterize w empirically, usually using 
two or more free parameters. Among all the parametrization forms of EOS, we will consider the Chevallier-Polarski-Linder 
(CPL) model (Chcvallie r and Pola rski 200 ll iLinder 2 003 ) . which is probably the most widely used parametrization, since it 
presents a smooth and bounded behavior for high redshifts, and a manageable two-dimensional parameter space. However we 
will also consider novel parametrizations recently introduced in ( |Ma fc Zhang 20l7| and ( |Lazkoz, Salzano and Sen dra 2010) 
to avoid the future divergency problem in the CPL parametrization, which turned out to be able to genuinely cover many 
scalar-field models, as well as other theoretical scenarios. It is worth noticing that all the dark energy models considered 
so far are in agreement with the observational data. As a consequence, unless higher precision probes of the expansion rate 
and the growth of structure are developed, these different approaches cannot be discriminated. This degeneration suggests 
a kinematical approach to the problem of cosmic acceleration, relying on quantities that are only weakly model dependent. 
The cosmographic approach is only related to the derivatives of the scale factor and it makes it possible to fit the data on 
the distance - redshift relation without any a priori assumption on the underlying cosmological model. It is based on the only 
assumption that the metric is spatially homogeneous and isotropic. The SNIa Hubble diagram extends up to z = 1.7 thus 
invoking the need for, at least, a fifth order Taylor expansion of the scale factor in order to give a reliable approximation of the 
distance - redshift relation. As a consequence, it could be, in principle, possible to estimate up to five cosmographic parameters, 
(h, 50, Jo, so, la) , although the still too small data set available does not allow to get a precise and realistic determination of 
all of them. Once these quantities have been determined, one could use them to put constraints on the dark energy models. 
We are reverting the usual approach that attempts to derive the cosmographic parameters as a sort of byproduct of the 
assumed theory. Here, we use the cosmographic parameters to parametrize the quantities that characterize the model so that 
each dark energy model is characterized by the same set of parameters (h, qo, jo, so, lo)- For constraining the cosmographic 
parameters, we use the Union2 Type la Supernovae (SNIa) data set, the Hubble diagram constructed from some Gamma Ray 
Bursts luminosity distance indicators, and gaussian priors on the distance from the Baryon Acoustic Oscillations (BAO), and 
the Hubble constant h (such priors have been included in order to help break the degeneracies among model parameters). 
Actually, observations of the type la supernovae are consistent with the assumption that the observed accelerated expansion is 
due to the non zero cosmological constant. However, so far the type la supernovae have been observed only at redshifts z < 2, 
while in order to test if w is changing with redshift it is necessary to use more distant objects. New possibilities opened up 
when the Gamma Ray Bursts have been discovered at higher redshifts, the present record is at z = 8.26 (|Greiner et al. 2009(1 . 
GRBs are however enigmatic objects. First of all the mechanism that is responsible for releasing the incredible amounts of 
energy that a typical GRB emits is not yet known (see for instance Meszaros 2006 for a recent review). It is also not yet 
definitely known if the energy is emitted isotropically or is beamed. Despite of these difficulties GRBs are promising objects 
that can be used to study the expansion rate of the universe at high redshifts ( Bradley 2003 ; Schacfcr 2003 ; Dai ct al. 2004 
IBloom et al. 20031 IFirmani at al. 20051 ISchaefer 20071 ILi et al. 20081 lAmati et al. 20081 ITsutsui et al. 2009)1 . Actually even if 
the huge dispersion (about four orders of magnitude) of the isotropic GRB energy makes them everything but standard candles, 
it has been recently empirically established that some of the directly observed parameters of GRBs are correlated with their 
important intrinsic parameters, like the luminosity or the total radiated energy, allowing to derive some correlations, which 
have been tested and used to standardize GRBs and to calibrate these relations, and to derive their luminosity or radiated 
energy from one or more observables, in order to construct a GRBs Hubble diagram. It has been shown that such procedure 
can be implemented without specifying the cosmological model, see for instance, (jpemianski, Picdipalumb o and Rubano 201 1[ 
Dcmianski & Picdipalumbo 2011 ) and references therein. In our analysis we use two GRB HD data sets: one sample consists of 
109 high redshift GRBs and has been constructed from the Amati E P] i - Ei ao correlation (here is the peak photon energy 
of the intrinsic spectrum and E{ BO the isotropic equivalent radiated energy), applying a local regression technique to estimate, 
in a model independent way, the distance modulus from the recently updated Union SNIa data set. The second GRBs HD 
sample is constructed from 66 Gamma Ray Bursts (GRBs) derived using only data from their X-ray afterglow light curve. To 
this end, we used the recently updated Lx - T a correlation between the break time T a and the X-ray luminosity Lx measured 
at T a calibrated (using SNIa) from a sample of Swift GRBs (Cardonc ct al . 2010|) . It is worth noting that such GRBs HD 
are based on the use of a single correlation and contain a statistically meaningful number of objects. The use of the Amati, 
E Pt i - Eisc, and the Lx - T a correlations then avoids the need of combining different correlations to increase the number of 
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GRBs with a known distance modulus. Since each correlation is affected by its own possible systematics and characterized by 
different intrinsic scatter so that combining all of them in a single HD can introduce unexpected features and hence bias the 
constraints on the cosmology. It turns out that these data sets are sufficient for our aim of testing and comparing the new 
parametrizations. Moreover in order to check if the results of our cosmographic analysis are biased due to the procedure used 
to calibrate the GRBs correlations we performed a consistency test: we actually apply a full bayesian approach, extracting, 
at the same time, the correlation coefficients and the cosmological parameters of the model from the observed quantities. 
Also to accomplish this task, we use the Markov Chain Monte Carlo simulations and compute, simultaneously, the full 
probability density functions (PDFs) of all the parameters of interest. This approach does not require any prior information 
on the cosmological model and yields results that are not plagued by any of the various limitations known in the literature, 
see for instance (Dcmianski, Picdipalumbo and Rubano 2011; Dcmianski & Picdipalumbo 2011). Since such a procedure is a 
demanding job from the point of view of computation time, it has been applied for the Amati relation only. The results turned 
out to be fully statistically consistent with the ones obtained by performing a local regression technique, thus indicating that 
it is not affected by any systematic bias induced by the calibration procedure. For the other data set we use the Markov Chain 
Monte Carlo simulations just to perform the cosmological tests. 

The scheme of the paper is as follows. In Section 2 we describe the basic elements of the cosmographic approach and 
explicitly derive series expansions of the scale factor and other relevant parameters. In Section 3 we describe the observational 
data sets that are used in our analysis. In Section 4 we describe some details of our statistical analysis and present results 
on cosmographic parameters obtained from three sets of data. In Section 5 we present constrains on dark energy models that 
can be derived from our analysis. General discussion of our results and conclusions are presented in Section 6. 



2 THE COSMOGRAPHY APPROACH 

Recently the cosmographic approach to cosmology gained increasing interest for catching as much information as possible 
directly from observations, retaining the minimal priors of isotropy and homogeneity and leaving aside other assumptions. 
Actually, the only ingredient taken into account a priori in this approach is the FLRW line element obtained from kinematical 
requirements 



ds 2 = —c 2 dt 2 + a 2 (t) 



1 — kr' 2 



(1) 



where a(t) is the scale factor and k = +1,0, —1 is the curvature parameter. Using this metric, it is possible to express the 
luminosity distance dc as a power series in the redshift parameter z, the coefficients of the expansion being functions of the 
scale factor a(t) and its higher order derivatives. This expansion leads to a distance - redshift relation which only relies on 
the assumption of the FLRW metric thus being fully model independent since it does not depend on the particular form of 
the solution of cosmic evolution equations. To this aim, it is convenient to introduce the following cosmographic functions 
IjVisser 2004|) : 
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When evaluated at the present time to these functions correspond to the cosmographic parameters, which are usually referred 
to as the Hubble, deceleration, jerk, snap and lerk parameters, respectively^]. Furthermore, it is possible to relate the derivative 
of the Hubble parameter to the other cosmographic parameters : 

H=-H 2 (l + q), (7) 

H=H 3 (j + 3q + 2) , (8) 

d 3 H/dt 3 =H 4 (s - 4j - 3q(q + 4) - 6) , (9) 

d 4 H/dt 4 = H 5 {I - 5s + W(q + 2)j + 30(g + 2)q + 24) , (10) 

1 Note that the use of the jerk parameter to discriminate between different models was also proposed in ( Salmi ct al. 2003} in the context 
of the statefinder parametrization. 
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where a dot denotes derivative with respect to the cosmic time t. With these definitions the series expansion to the 5th order 
in time of the scale factor is: 



■^1=1 + H (t - to) - ^H 2 (t - t ) 2 + 3 ±H 3 (t - to) 3 + S ^H A o{t - t f + l ^H 5 (t - t ) 5 + 0[(t - t f] . 



(11) 



D = c / dt = c(t - t») , 

we can construct the series for z(D), actually 
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To obtain the physical distance D expressed as a function of redshift z we reverse the series z(D) — > D(z), obtaining: 



D{z) = ^- (V° + Viz + Viz 2 + V\z A + Vjz 4 + 0(z 5 )) 
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In the following we will be interested not in the physical distance D(z), but in the luminosity or angular-diameter distance, 
which can be calculated as 
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where ro(D) is: 
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sinh(f to D £-M)k = -1. 
If we insert the series expansion of a(t) in ro(D), we have the cosmographic expansion of ro(D): 
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In our analysis we will consider spatially flat cosmological models only, so that 
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With this expansion the luminosity distance is given as: 
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While for the angular diameter distance we get: 
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It is worth noting that one can obtain the same final expression for the distance starting from the Taylor series expansion of 
the Hubble parameter instead of the scale factor, namely: 



H( z y 
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To compute all the terms we use the derivation rule 
d 



-(l + z)H- 



(49) 



d_ 

dt y ' ' dz 

The series expansion of the Hubble parameter will be used in our analysis in the definition of our Markov chains algorithm given 
the observational data we use. It is worth noting that since the cosmography is based on series expansions, the fundamental 
difficulties of applying such an approach to fit the luminosity distance data using high redshift distance indicators are connected 
with the convergence and the truncation of the series. Recently the possibility of attenuating the convergence problem has 
been analyzed by defining a new redshift variable, see ( |Vitagliano et al. 2010[ ), the so called y-redshift: 
z 



*^y = — z . (50) 

It turns out that for a series expansion in the classical z-redshift the convergence radius is equal to 1, which is a drawback 
when one wants to extend the application of cosmography to redshifts z > 1. The y-redshift could potentially weaken this 
problem because the z-interval [0, oo] corresponds to the y- interval [0, 1], so that we are mainly inside the convergence interval 
of the series, even for CMB data (z = 1089 — > y = 0.999). Thus, in principle, we could extend the series up to the redshift of 
decoupling, and one could place CMB related constraints within the cosmographic approach 0. However even using the series 
expansions in y-redshift the problem of the truncation of the series remains. Here we consider a forth order expansion and 
are able to successfully put bounds on the parameters in a statistically consistent way. In order to give reasonably narrow 
statistical constraints we apply a Markov Chain Monte Carlo (MCMC) method, which allows us to obtain marginalized 
likelihoods on the series coefficients from which we infer rather tight constraints on those parameters. Actually, in our code we 
have inserted several tests, which give us control over several physical requirements we expect from the theory. For instance, 
since we use data related to the Hubble parameter H(z), we are able to set restrictions on the Hubble parameter, Ho = H(0), 
and thus to obtain a considerable improvement in the quality of constraints. 



3 OBSERVATIONAL DATA SETS 

In our cosmographic approach we use the currently available observational data sets on SNIa and GRB Hubble Diagrams, 
and we set gaussian priors on the distance from the Baryon Acoustic Oscillations (BAO), and the Hubble constant h. Such 
priors have been included in order to help break the degeneracies among the parameters of the cosmographic series expansion 
in Eqs. (p6) . 



3.1 Supernovae 

Over the last decade the confidence in type la supernovae as standard candles has been steadily growing. Actually, the SNIa 
observations gave the first strong indication of an accelerating expansion of the universe, which can be explained by assuming 
the existence of some kind of dark energy or nonzero cosmological constant. Since 1995 two teams of astronomers - the High-Z 
Supernova Search Team and the Supernova Cosmology Project - have been discovering type la supernovae at high redshifts. 
First results of both teams were published by (Rie ss et al. 1998)) and (fPerlmutte r et al. 1999[) . Here we consider the recently 
updated Supernovae Cosmology Project Union2 compilation (Amanullah et al. 2010), which is an update of the original 
Union compilation, now bringing together data on 719 SN, drawn from 17 data sets. Of these, 557 SN, spanning the redshift 
range (0.015 < z < 1.55), pass usability cuts and outliers removal, and form the final sample used to constrain our model. 
We actually compare the theoretically predicted distance modulus /i(z) with the observed one, through a Bayesian approach, 
based on the definition of the distance modulus, 

fi(z 3 ) = 51og 10 (£> L fe, {#,})) + fi , (51) 

where Dl{zj, {9i}) is the Hubble free luminosity distance, expressed as a series depending on the cosmographic parameters, 
0i = (qo,jo, so, lo)- The best fits were obtained by minimizing the quantity 



2 Let us note that the introduction of this new redshift variable will not affect the definition of cosmographic parameters. 
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3 = 1 a ^ 



(52) 



where the j are the measurement variances. The parameter /j,o encodes the Hubble constant and the absolute magnitude 
M, and has to be marginalized over. Giving the heterogeneous origin of the Union data set, and the procedures for reducing 
data, we have worked with an alternative version of Eq. I|52[). which consists in minimizing the quantity 

Xsn(W) =<*--, (53) 

C3 

with respect to the other parameters, where 
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It turns out that xIjv ls j us t a version of Xsjvi minimized with respect to /io- Actually, we find that 

Xsn(mo, {9i}) = ci - 2c 2 ^ + c 3 Mo , (57) 

which clearly becomes minimum for fio = C2/C3, and so we see that xIn = xInCa'o = C2/C3, Furthermore, one can check 

that the difference between Xsn an< l Xsn is negligible. 



3.2 GRBs Hubble diagram 

By virtue of their enormous energy release, GRBs are visible up to very high z, and hence are ideal candidates for our high- 
redshift cosmography task. Unfortunately, GRBs are everything but standard candles because their peak luminosity spans a 
wide range. Nevertheless there have been many attempts to make them standardizeable candles resorting to the use of some 
empirical correlations among distance dependent quantities and rest frame observables ( |Amati et al. 2008]) . Such empirical 
relations allow one to infer the GRB rest frame luminosity or energy from an observer frame measured quantity so that 
the distance modulus can be estimated with an error mainly depending on the intrinsic scatter of the adopted correlation. 
Combining the estimates from different correlations, Schaefer (2007) first derived the GRBs HD for 69 objects, which has 
been further enlarged using updated samples, different calibration methods and also different correlation relations, see for in- 
stance (Cardonc et al. 2008), (Dcmianski, Picdipalumbo and Rubano 2011), (Dcmianski & Picdipalumbo 2011), showing the 
interest in the cosmological applications of GRBs. In this paper we perform our cosmographic analysis using two GRBs HD 
data set, build up by calibrating the Amati E p ,\ - E{ so and the Lx - T a correlations respectively. 



3.2.1 The calibrated Amati Gamma Ray Bursts Hubble diagram 

It has been recently empirically established that some of the directly observed parameters of Gamma Ray Bursts are connected 
with the isotropic absolute luminosity Li SO , the collimation corrected energy E 1 , or the isotropic bolometric energy i?i so of a 
GRB. These quantities appear to correlate with the GRB isotropic luminosity, its total collimation-corrected or its isotropic 
energy. The isotropic luminosity cannot be measured directly but rather it can be obtained through the knowledge of either 
the bolometric peak flux, denoted by Pboio, or the bolometric fluence, denoted by Sboio- The isotropic luminosity is given by 

L iso = 47rdl(z)P bolo , (58) 

the total collimation-corrected energy reads as 

E iao =47Tdl(2)Sbolo(l + 2)" 1 , (59) 

and the total collimation-corrected energy is 

E-, = 47rdl(2)SboioF bcam (l + z)- 1 , (60) 

where -Fbcam is the beaming factor. 

Therefore, Li SO , E^ and -Bi so depend not only on the GRB observables Pboio or SWo, but also on the cosmological 
parameters, through the luminosity distance di(z). As a consequence, there is a big circularity problem to overcome, since it 
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Figure 1. Distance modulus /J.(z) for the calibrated GRBs Hubble diagram made up by fitting the Amati correlation. 



is not immediately possible to calibrate such GRBs empirical laws, and to build up a new GRBs Hubble diagram, without 
assuming any a priori cosmological model. In ( [Demianski fc Pi cdipalumbo 2 011[ ) we have applied a local regression technique 
to estimate, in a model independent way, the distance modulus from the recently updated Union SNIa sample, containing 
557 SNIa spanning the redshift range of 0.015 < z < 1.55. The derived calibration parameters have been used to construct 
an updated GRBs Hubble diagram. In particular, by using such a technique, we have fitted the so-called Amati relation and 
constructed an updated Gamma Ray Bursts Hubble diagram, which we call the calibrated GRBs HD, consisting of a sample 
of 109 objects, shown in Fig. [1] Their redshift distribution covers a broad range of z, from 0.033 to 8.26, thus extending 
far beyond that of SNIa (z max ~ 1-7), and including GRB 092304, the new high-z record holder of Gamma Ray Bursts. 
Since here we want to use such calibrated GRBs HD to perform the high-redshift cosmographic analysis described below, it 
is worth to discuss some arguments about the reliability of using the Amati relation for cosmological tasks. For instance, one 
of the discussed objection is a supposition that the Amati relation is dominated by various selection effects of the detectors 
and differences within the GRB population. Recently such a thesis has been developed according to a conceptually simple 
argument based on the fact that the Si, ;o--E P cak diagram presents two limit lines, where bursts cannot be below if the Amati 
and Ghirlanda relation holds (Nakar & Piran 2005), (Collazzi et al. 2012"]) ■ Actually it turns out that 

£ = — =#"10" — -, (61) 

oboio (i _|_ z y + ~ 

where K and rj are the normalization and the slope of the Amati corralation (E pea k = KE^ so ), and a corresponds to 
the scatter of the data points around the rest frame correlation. Let us note that the left side of the Eq. Q6ip uses only 
directly observable quantities, while the right side is only a function of distance. As the distance rises, d 2 L gets larger and 
(1 + z)~( 1+v ' gets smaller, giving rise to a maximum value for the right side, which cannot exceed therefore a limit value. 
Indeed an immediate test of the Amati relation can be performed, by investigating where the average burst falls below the 
Amati limit. In (ICollazzi et al. 2005|) the results of such an analysis have been interpreted as a clear disproof of the Amati 
relation. However such a conclusion sounds premature, and should be postponed untill the available datasets will be more 
consistent and once some sources of uncertainty have been taken under control. For instance, the limit lines could vary with 
the background cosmology, or with the slope parameter, r/, in a not negligible way, as indicated by Fig. [2] and Fig. [3] with 
sensitive effects on the fraction of violators of such a limit ( on which the analysis itself is based). 

Here we do not intend to analyse thoroughly such aspects, which we postpone to a forthcoming paper, and which, even if 
important, do not touch the heart of our cosmographic investigation; however we want just to note, as an a posterior argument, 
that if we relay on the Amati relation, and use it to build up the Hubble diagram, it turns out that it is fully consistent with 
the HD obtained from all the other relations available for the GRBs, as seen in Figs. [4] In Fig. [5] we show how such limit 
lines can be shifted, by varying the energy ratio as in the Figure [2] It turns out that the number of outliers strongly depends 
on the position of the limit lines. Moreover, in (Nava & al. 2 011jl the spectral properties of short and long GRBs, detected 
by the Gamma-ray Burst Monitor ( GBM), have been studied over an unprecedented wide energy range. It turns out that the 
fraction of long GRBs, which are outliers (at more than 3cr) with respect to the Amati relation, is ~ 3 per cent, while there 
are no outliers (at more than 3a) for the E pea k — Li SO correlation. 
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y 



Figure 2. Predicted value of the Amati relation energy ratio £ as a function of redshift, varying the cosmological parameters and the 
slope of the Ei ao -E pe ^. Here is illustrated an example of the ACDM model where the Hubble constant h varies from h = 0.74 to h = 0.65 
(narrow blue region), and the slope rj _1 from n^ 1 = 1.8 to n^ 1 = 2.18 (large pink region). It turns out that the peak can be shifted 
noticeably. 




to- 8 
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Figure 3. The basics of the Nakar and Piran test in a graphical form. Any burst (even without a known redshift) can be plotted on 
this diagram. If the Amati relation is correct, then all burst should lie above the solid line. Here it is shown how such a limit value can 
be shifted, by varying the energy ratio as in the previous Figure. The superposed points are taken from Wei 2010 ; it turns out that the 
number of outliers strongly depends on the position of the limit lines. 

3.2.2 The Lx - T a Gamma Ray Bursts Hubble diagram 

The Lx - T a correlation between the luminosity Lx at the break time T a and T a itself is the only empirical law relating 
quantities measured from the afterglow light curve, as described by the universal fitting function proposed by Willingale et al. 
(2007, hereafter W07), rather than being related to the prompt emission quantities. It has been first discovered by Dainotti 
et al. (2008) and later confirmed by the semiempirical models of Ghisellini et al. (2009) and Yamazaki (2009). More recently, 
Dainotti et al. 2010 ( herefater D10) have increased the GRBs sample and rederived the Lx-T a correlation, selecting a class 
of high luminosity long GRBs with very well measured (Lx,T a ) parameters and lightcurve closely matching the W07 model. 
Referring to this class of objects as canonical GRBs, D10 have demonstrated that they define an upper envelope for the 
Lx-Ta correlation with the same slope, but a higher intercept than that for the full sample. Cardone et al. (2011) used 
this Lx-Ta Gamma Ray Bursts Hubble diagram to constrain cosmological parameters of some simple dark energy models. 



40 42 44 46 48 40 42 44 46 48 

A'Amati Another 

Figure 4. Left panel: comparison between the HD obtained with the Amati relation and the L — ri ag , L — V , L — E pea k, L — trt 
relation. Right panel: behaviour of the residuals among the Amati and the other relations. It turns out that these datasets are fully 
consistent and strongly correlated with the Spearmans correlation p = 0.8. 
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Figure 5. Left panel: comparison between the HD obtained with the Amati relation and the Lx -T a relation. Right panel: behaviour 
of the residuals. It turns out that these data sets are fully consistent and strongly correlated with the Spearmans correlation p = 0.72. 




Figure 6. Distance modulus fi(z) for the calibrated GRBs Hubble diagram made up by fitting the Lx — Ta correlation. 



Since the canonical sample is statistically poor (it consists of only 8 objects), here we use the full sample to perform our 
cosmographic analysis. We build up such sample calibrating, the u < 4 sample from D10, with Local Regression technique 
used in (Dcmianski, Picdipalumbo and Rubano 2011) and (Dcmianski & Picdipalumbo 2011), for the 



log L x =a log ( yq~) + b 



+ «4 



correlation. Here u is an error parameter: u — ^Jo\ 
simply note that Lx is related to the luminosity distance oIl(z) as 

Lx=4*dl( Z ){l + z)- (2+p) F x , 



(62) 



In order to infer the distance modulus of each GRB, we then 



(63) 



ft being the slope of the energy spectrum (modelled as a simple power - law) and Fx the observed flux both measured at the 
break time T a . Having measured (T a , f3, Fx) and inferred Lx using Eg . (f52"|) . we can then estimate the GRB distance modulus 
as : 



fi(z) = 25 + 51ogd L (z) = 25 



■log 



4tt(1 + z)-( 2 +«Fx 



(64) 



where oIl(z) is in Mpc. Such (Lx, T a )-GRBs Hubble diagram in shown in Fig. [6] The uncertainty is estimated by propagating 
the errors of (/3, Fx , Lx)- While both SNIa and GRBs are based on the concept of standard candles, an alternative way 
to probe the background evolution of the universe relies on the use of standard rulers. Nowadays the Baryonic Acoustic 
Oscillations (BAOs) which are related to the imprint of the primordial acoustic waves on the galaxy power spectrum are 
widely used as such rulers. In order to use BAOs as constraints, we follow (Percival a l. 2010[) by first defining: 

with Zd the drag redshift computed using the approximated formula in (Eiscnstein and Hu 1998), r s (z) the comoving sound 
horizon given by : 



a 2 H(a)y/l + (3/4)fi 6 /fi 7 
and Dv(z) the volume distance defined by (Eiscnste in et al. 2 005) : 
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D v (z) = 



H[z) 



Dl{z) 



1 + 2 



1/3 



(67) 



4 STATISTICAL ANALYSIS 



In this section we describe some details of our statistical analysis and present our main results on the constraints for the 
cosmographic expansion parameters from the current observational data sets described above. In order to constrain the 
cosmographic parameters, we perform a preliminary and standard fitting procedure to maximize the likelihood function 
£(p) oc exp [— x 2 (p)/2], where p is the set of cosmographic parameters and the expression for X 2 (p) depends on the data set 
used. As a first test we consider only the SNIa data, thus we define : 

E 



x 2 (p 



Vobs{Zi) - fJ, t h(Zi, p) 



h - 



n.742y 



+ 



0. 1.356 y 



(68) 



Here. 



0.036 j \ 0.0034 

fiobs and nth are the observed and theoretically predicted values of the distance modulus, while the sum is over all 



the SNIa in the sample. The last two terms are Gaussian priors on h and um = QmIt 2 and are included in order to help 
break the degeneracies among the model parameters. To this aim, we have resorted to the results of the SHOES collaboration 
(R.ics s~et al. 2009]) and the WMAP7 data l|Komatsu et al. 2010 jl . respectively, to set the numbers used in Eqs. (|68p . When we 
are using GRBs only, we define : 

Ngrbhd 



x 2 (p)= E 



flobs(Zi) — Hth(zi, p) 



( h- 0.742 \" n 
\ 0.036 ) + V 



* s 1 356 ^ 2 



£(p)oc 



0.036 / V 0.0034 
As a next step, we combine the SNIa and GRBs HDs with other data redefining £(p) as : 

CX P (~XsNIa/GRB/2) 



(69) 



(27T) 



SNIa/ G F! B 



exp 



Csnici/grbI 1 / 
1 ( h - h ob . N : 

Oh 



exp(-xlAo/2) 



(27r)^Ao/2| CsAO |i/2 



exp 



exp(- X i/2) 



n-n ob 

OR 



(70) 



( 27r )Af„/2j Cff |l/2 ' 

The first two terms are the same as above with CsNia/GRB the SNIa/GRBs diagonal covariance matrix and (h t s ,Oh) = 
(0.742,0.036). The third term takes into account the constraints on d z — r a (zd)/Dv{z) with r a (zd) the comoving sound 
horizon at the drag redshift Zd (which we fix to be r a (zd) = 152.6 Mpc from WMAP7) and the volume distance is defined as 
in Eq. (|67p . The values of d z at z — 0.20 and z = 0.35 have been estimated by Percival et al. (2010) using the SDSS DR7 
galaxy sample so that we define Xbao ~ I-^C^qC with D T = (c^ 6 ! — d^, do'ls — dQ% 5 ) and Cbao is the BAO covariance 
matrix. The next term refers to the shift parameter HBond et al. 199 7 Efstathiou & B ond 1999)) : 

dz' 



11 = # A 



H(z') ' 



(71) 



with Zi, = 1090.10 the redshift of the last scattering surface. We follow again the WMAP7 data setting (1Z b s ,on) = 
(1.725, 0.019). While all these quantities (except for the Gaussian prior on h) mainly involve the integrated E(z), the last term 
refers to the actual measurements of H(z) from the differential age of passively evolving elliptical galaxies. We then use the 
data collected by Stern et al. (2010) giving the values of the Hubble parameter for Mh = 11 different points over the redshift 
range 0.10 < z < 1.75 with a diagonal covariance matrix. We finally perform our cosmographic analysis, considering a whole 
data set containing both the SNIa Union data set and the calibrated GBRs HD (which we call the cosmographic dataset), and 
slightly modifying the likelihood £(p). To compute the likelihood, we use a Bayesian approach, implementing a Monte Carlo 
Markov Chain technique. For each Monte Carlo Markov Chain calculation, we run two times four independent chains that 
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Parameter 


h 


90 


jo 


so 


lo 


Best Fit 


0.755 


-0.462 


0.134 


0.696 


-3.57 


Mean 


0.761 


-0.471 


0.258 


0.164 


-0.453 


2 cr 


(0.71,0.81) 


(-0.53, -0.40) 


(0.105,0.68) 


(-0.996, 1.28) 


(-8.23,7.37) 



Table 1. Constraints on the parameters of Cosmography I (from combining the SNIa HDs with BAO and H(z) data set (2<r error bars)). 



Parameter 


h 


?o 


30 


so 


lo 


Best Fit 


0.742 


-0.406 


0.86 


29.767 


-7.757 


Mean 


0.732 


-0.359 


0.28 


20.33 


1.044 


2 cr 


(0.67,0.81) 


(-0.53, -0.134) 


(0.104,0.95) 


(-0.86,29.72) 


(-30.01,29.95) 



Table 2. Constraints on the parameters of the Cosmography II (from combining the SNIa HD, the Amati Gamma Ray Bursts HD with 
BAO and H(z) data sets (2cr error bars)). 



consist of about 100000 chain elements each q We test the convergence of the chains by the Gelman and Rubin criterion, 
finding R — 1 of order 0.01, which is more restrictive than the often used and recommended value R — 1 < 0.1 for standard 
cosmological investigations. Moreover in order to reduce the uncertainties on cosmographic parameters, since methods like 
the MCMC are based on an algorithm that moves randomly in the parameter space, we a priori imposed some constraints 
on the series expansions of H 2 (z) and cIl(z). The most general and obvious constraint is the positivity requirement^: 

• d L (z) > , 

• H 2 {z) > , 

applied for all our redshift ranges. We first run our chains to compute the likelihood in Eqs. (|68[l and/or (I69|l . using as starting 
points the best fit values obtained in our pre- statistical analysis, in order to select more efficiently the space (cosmographic) 
parameters region, and mainly to select the starting points. Therefore we perform the same Monte Carlo Markov Chain 
calculation to evaluate the likelihood in Eq. (|70[) , combining the SNIa HD, the BAO and H(z) data with the GRBs HD 
respectively, as described above. We will refer to such two different cases as Cosmography I and Cosmography II. we throw 
away first 30% of the points iterations at the beginning of any MCMC run, and we thin the two-runned chains. We finally 
extract the constraints on the parameters, coadding the thinned chains. Both Cosmography I and Cosmography II are 
implemented using z (z-Cosmography I /II) and y (y-Cosmography I/II) series. In Table [T] we present the results of our 
cosmographic analysis applied to the Union2 SNIa data set. It turns out that in Cosmography I the snap parameter, so, 
is weakly constrained and the lerk, lo, is actually unconstrained. It turns out that with the y-Cosmography I we obtain 
practically the same results, without any statistical meaningfull difference. However the j/-redshift approach allows a faster 
convergence of the chains. 

In the same way, the z-Cosmography II and y-Cosmography II turned out to be fully statistically equivalent and in Table 
[2]we present the results of our cosmographic analysis adding the calibrated Amati Gamma Ray Bursts Hubble diagram in the 
y-redshift case (y-Cosmography II). In Fig. [7] we plot the marginalized likelihood function for the deceleration parameter qo- It 
turns out that in this Cosmography II the snap parameter, so, is weakly constrained and the lerk, lo, is actually unconstrained. 

It is worth noting that from our statistical MCMC analysis it turns out that the deceleration parameters qo is clearly 
negative in all the cases. The marginal likelihood distribution for the current deceleration parameter qo indicates that there 
is just a negligible probability for qo > 0, as shown in Fig. [Jj Moreover the value of the jerk jo is significantly different from 
the ACDM value jo ~ 1, what is also indicated by the marginal likelihood distribution, see Fig. [8] In Fig. [9] are shown the 
confidence regions for h, qo, and jo- the left-side and the right-side panels concern the (h-qo) and the (go-jo) plane, respectively. 
In Fig. [10] we plot the observational data compared with the maximum likelihood curve. 

Independently of the cosmic deceleration today, it is of interest to investigate if there is another change in the sign of the 
cosmic acceleration after a prior transition from a decelerated to an accelerated phase at moderate redshifts (z t ~ 0.5 — 1). We 



3 Actually multiple independent chains can be started in different points of the parameter space to ensure good mixing, i.e. an adequate 
exploration of the whole parameter space. 

4 Let us note that here we do not use the constraint < Q m < 1 , that was employed, for instance, in 
( Capozzicllo, Lazkoz, and Salzano 2011), since it is a sort of meta cosmographic constraint, since it would require to postulate the 
form of the Friedman equations. However, when in the following section we are going to investigate the implications of our results on 
the evolution of the dark energy equation of state for some parametrizations, we use such a constraint as an a posteriori control. 
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Figure 7. The marginalized probability density function (PDF) for the deceleration parameter qo, as provided by Cosmography II. 
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Figure 8. The marginalized probability density function (PDF) for the jerk parameter jo, as provided by Cosmography II. 




Figure 9. Confidence regions in the (h-qo) and the (qo-jo) plane, respectively, as provided by Cosmography II. The inner brown region 
define the 3<x confidence level. It turns out that the parameters h, qo and jo are well-constrained, that the values qo > are ruled out, 
the value jo = 1 (which is the ACDM value) is statistically not favourable. 
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Figure 10. Comparison of the observational data used in our Cosmography II with the maximum likelihood curve (the thick red curve) 
visible in the center of the filled region, corresponding to the 2c confidence levels for the cosmographic parameters. 
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Figure 11. Reconstruction of deceleration history. The filled area delimits the 2cr confidence region for the q(z) reconstruction obtained 
from our Cosmography II. 
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Figure 12. Reconstruction of deceleration history obtained from our Cosmography I. 



use the cosmographic parameters (go, jo, so and Zo) to reconstruct q(y) and to trace the deceleration history of the universe. 
Actually, from the power series expansion of the scale factor one can also express the deceleration parameter as a power series 
in time, according to the definition in Eq. [2] This time dependent parameter can be written as a power series in y-redshift, 
using the derivation rule 

Tt=^^ H Ty- (72) 

The derivation of a power expansion (of fourth order) for q(y) from the scale factor expansion allows for a decelerated past, a 
transition to an accelerated phase, a point of maximum acceleration, then a slowing down of the acceleration and a transition 
to a recent or future decelerating phase, as shown in Fig. 1111 Appearance of transient acceleration is predicted or allowed 
by several dynamic models ( Shaficloo, Sahni, and Starobinsky 2009). In contrast, the ACDM model predicts a monotonic 
deceleration history connecting its asymptotic limits in the past and future, q(z — > oo) = 0.5 and q(z — » —1) = — 1. The 
q(z) reconstruction obtained from our cosmographic parametrization is shown in figure 1111 It is worth noting that also the 
q(z) reconstruction obtained from our Cosmography I allows a transient acceleration, as shown in Fig. 1121 but the specific 
properties of the deceleration history of the universe is rather different from the one obtained if we include GRBs HD in our 
cosmographic parametrization. 



4.1 Cosmography with the Lx - T a Gamma Ray Bursts Hubble diagram 

In this section we shortly describe the results of our statistical analysis of another data set formed by combining the BAO and 
H(z) data with the Lx-T a Gamma Ray Bursts Hubble diagram described above (referred in the following as Cosmography 
III). We have decided not to include such GRBs data set in the overall analysis, performed in Cosmography II, by virtue of 
the uncertain features of this recently discovered correlation, which should be further investigated in order to get stronger 
confidence and confirmation. It turns out that the results are mostly compatible with the previous ones obtained in our 
Cosmography I and II, as indicated in Tab. [3] Moreover it is worth noting that also in this case the value jo = 1 (which is 
the ACDM value for the jerk) is statistically not favourable, as shown also in Fig. 1131 

The q(z) reconstruction allows a transient acceleration, shown in Fig. 1141 as in Cosmography I and II, strengthening the 
reliability of the Lx - T a correlation. 
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Parameter 


h 


90 


jo 


so 


l 


Best Fit 


0.74 


-0.44 


0.386 


-0.719 


3.82 


Mean 


0.713 


-0.402 


0.45 


2.54 


-24.1 


2 a 


(0.68,0.79) 


(-0.50, -0.26) 


(-0.73, 1.4) 


(-1.3,8.4) 


(-72.0,8.5) 



Table 3. Constraints on the parameters of the Cosmography III (from combining the SNIa HD, the Lx-T a Gamma Ray Bursts HD 
with BAO and H(z) data sets (2<r error bars)). 
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Figure 13. Confidence regions in the (h-qo) and the (qo-jo) plane, respectively, as provided by the Cosmography III. It turns out that 
these results are compatible with the previous ones obtained in our Cosmography II. 



5 IMPLICATIONS FOR DARK ENERGY 

In this section we investigate the implications of the results from our simulations on cosmography for different and specific 
parametrized dark energy models. The link between the cosmographic and the dark energy parametrization is based on the 
series expansion (in redshift) of the Hubble function H(z). Actually, for a spatially flat cosmological model it turns out that: 



H{z) = H oy /{\ - n m )g(z) + Qm(z + l) 3 , 
H d (z) = -(z + l)H(z)H'{z), 
H 2d (z) = -(l + z)H(z)H' d (z) y 
H 3d (z) = -(l + z)H(z)H' 2d (z), 
H 4d (z) = -(l + z)H(z)HUz), 

3 f ■"(*) + ! dx 

where g(z) — exp Jo x + 1 ; and w(z) any parametrized form of the dark energy equation of state. It turns out that 
lim H d (z)=-H (l + qo) , 

z — >0 

lim H 2d (z)=Hl (jo + 3q + 2) , 

z->0 

lim H 3d (z)=Ho (s - 4j - 3q (q + 4) - 6) , 

z->0 

lim H u {z)=Ht (Jo - 5s + 10(g + 2)j + 30(g + 2)*, + 24) . 

z->0 

In this section we consider three different parametrizations: 
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Figure 14. Reconstruction of deceleration history obtained from our Cosmography III. The filled area delimits the 2cr confidence region 
for the q(z) reconstruction obtained from our analysis. 
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• the CPL parametrization for dark energy given by 

111(2) = wo + Wiz(l + zy 1 , (82) 

• a novel parametrization recently introduced in ( |Ma fc Zh ang 201l| to avoid the future divergency problem of the CPL 
parametrization, and to probe the dynamics of dark energy not only in the past evolution but also in the future evolution, 

, n / sin(z + 1) , A , 

w(z) =w + Wl I z + 1 -sin(l) J , (83) 

• an oscillating dark energy equation of state recently discussed in (?) 

w 1 (l-co S (8lo g (z + l))) +wo 
\og(z + 1) 



It is worth noting that the efficiency of such investigation, i.e the possibility of inverting the equations (|78H8ip . strongly depends 
on the number of cosmographic parameters we are working with and on how many DE equation of state parameters we are 
going to consider as free. For instance, restricting to the CPL equation of state given in Eq. 1821 one has three possibilities: 

• with two cosmographic parameters, (qo,jo), we can derive some information about a constant dark energy model (i.e. 
wi =0), with: 

2{j -q -2ql) 



fim(<70, jo) 



1 + 2 jo 



, . . 1 + 2j - 6g , OE , 

wo(«7o,Jo) = 5—7 , (85) 

—3 + ogo 

• with two cosmographic parameters, (qo,jo), we can derive some information also about a dynamical dark energy model, 
(i.e. Wi 7^ 0), leaving Q m free, with: 

l-2q 



wo(qo, fi m ) = 



3(-l + fi m ) 



wi(go,jo,n m )= - — (-2jo(-l + n m )-2g x (l + 2q ) + n m (-l + 6qo)) , (86) 

3(— 1 + \i m y 

• with three cosmographic parameters, (go, jo, so), we can derive some information about a dynamical dark energy model, 
with Q. m depending on the cosmographic parameters, i.e. Q m = Q m (qo, jo, so). The same holds true for wo = wo(qo, jo, so) 
and mi = w 1 (q ,jo,s ). 

For these relations all the statistical properties of these parameters (median, error bars, etc.) can be directly extracted from 
the cosmographic samples we have obtained from the MCMC analysis. In our investigation we prefer to adopt a conservative 
approach, considering only the parameters which are well constrained by our cosmographic analysis, that is (go, and jo) 
leaving Q m free. It is therefore possible to obtain wo and w\ in terms of fi m and other cosmographic parameters. For the 
parametrizations of the equation of states examined above, we get 

• parametrization given by Eq. (J83J) 

, . _ . 2j (ftm - 1) + 2g (-3fi m + 2g + 1) + £l„ 
u>i(g ,jo,fim) = 



3(fi m - l) 2 (cos(l) -sin(l)) 
parametrization given by Eq. (|84ll 



Wo(go '" m)= 3(n ro - i) (89) 

, os -2jo (Qm - 1) + (6go - 1) n m - 2q (2g + 1) 

m(q ,jo,n m )= w(n m -iy • (90) 

In Eq. [90] the parameter S can be expressed as a function of fi m and lo, inverting Eq. 1811 which in our case gives 

lo = - (3 (S 4 wi (On - 1) (3wi (n m - 7) + 2) + S 2 Wl (9w (fi m - 1) (w a (110 m - 23) + 8fi m - 31) - 92fi m ) + 

wo (fi m - 1) (3w (6 (6w (wo + 2) + 7)Q m - 3w {18w + 47) - 134) - 163)) + 2765 2 wi + 70) . (91) 

Actually, in Figs. 1151 1161 and [17] we can reconstruct the redshift behaviour of the different dark energy equations of state 
corresponding to different values of Q m (left panel), and obtained when the cosmographic parameters are varying within the 
2a region of confidence, and Q m is fixed (and set to be f2 m = 0.22). The cosmographic analysis allows us to infer the actual 
value of the parameters appearing in the dark energy equation of state, and, as expected, can provide constraints on its 
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Figure 15. Upper Panel The y-redshift dependence of the equation of state described by the CPL paramctrization (see Eq. I82| l for 
different values of the Q m , as described in the label. Bottom Panel The j/-redshift dependence of the CPL equation of state for different 
values of the cosmographic parameters qo and jo- The filled region corresponds to the allowed behaviour of the equation of state, when 
the cosmographic parameters are varying within the 2<r region of confidence. Q m is fixed and set to be f! m = 0.22. 



redshift evolution mainly in a low redshift range. It turns out that the equation of state is evolving for all the parametrization 
considered, as confirmed for example in Figs. 1181 where it is clear that for the parameter wi, appearing in the non constant 
term of the equation of state, the case wi = has marginal confidence in the 2a region of confidence for the parameters 
go, jo, independently of the value of Q m - Moreover, Fig. [19] reflects the possibility of a deviation from the ACDM cosmological 
model. Results of our cosmographic analysis are only marginally compatible with predictions of the ACDM model. In a 
forthcoming paper we are going to reconstruct the EOS of dark energy from the observational data. As a final remark we note 
that the reliability of our cosmographic analysis is strongly related to the question whether there is a relation between the 
highest expansion order in the Taylor series and the redshift range where this series can be applied. One could expect that 
we would need a series expansion truncated at higher orders when increasing up the redshift range. However, the errors on 
the cosmographic coefficients will increase when higher order expansions are considered. As those parameters are correlated 
among them, errors in the low order series coefficients propagate to the additional coefficients included in the higher order 
series. If these errors turn out to be too large, the strength of cosmography will vanish. In order to estimate the error resulting 
from stopping the expansion at the fifth order we analyze the (relative) residuals between the exact distance modulus and 
fifth order series expansions of the same quantity, when varying the cosmographic parameters go, jo, so, and lo within the 3a 
region of confidence, and for the CPL parametrization considered above, as illustrated in Fig. 1201 It turns out that actually 
the error can be significant (with respect to the data) already for y > 0.8 (corresponding to z ~ 4). As far as the results of 
our analysis are concerned, they are not dramatically affected by such truncation error, not only because few data-points falls 
into the forbidden region, but mainly since we expect to estimate at most go, jo said so (we actually consider lo essentially 
unbounded), that is the cosmography parameters connected to the expansion at forth order. Moreover when we investigate 
the implications of the results from our simulations on cosmography for different and specific parametrized dark energy model, 
we have considered only the parameters which are well constrained by our cosmographic analysis, that is go, and jo- Finally 
it is worth noting that Fig. [20] suggests that, because of the actual precision of the observations, cosmography should not be 
used to calibrate the GRBs correlation relation, even if one limits the procedure to GRBs at y < 0.6. 
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Figure 16. Upper Panel The j/-redshift dependence of the equation of state described by the parametrization (see Eq. [83} for different 
values of the fi m , as described in the label. Bottom Panel The y-redshift dependence of the same equation of state for different values 
of the cosmographic parameters go and The filled region corresponds to the allowed behaviour of the equation of state, when the 
cosmographic parameters are varying within the 2<r region of confidence. O m is fixed and set to be fi m = 0.22. 



6 DISCUSSION AND CONCLUSIONS 

In this paper we are studying the possibility to extract model independent information about the dynamics of the universe 
by using a cosmographic approach considering only minimal assumptions (isotropy, homogeneity, Taylor series expansion of 
distances) without choosing any dynamical model a priori. In order to explore it systematically, we performed an high-redshift 
analysis that allowed us to put constraints on the cosmographic expansion up to the fifth order, based on the Union2 Type la 
Supernovae (SNIa) data set, the Hubble diagram constructed from some Gamma Ray Bursts luminosity distance indicators, 
and gaussian priors on the distance from the Baryon Acoustic Oscillations (BAO), and the Hubble constant h. Actually we 
use two GRB HD data set: one sample consists of 109 high redshift GRBs and has been constructed from the Amati E P: i - 
Eiso correlation. The second GRBs HD sample is constructed from 66 Gamma Ray Bursts (GRBs) derived using only data 
from their X-ray afterglow light curve. To this end, we used the recently updated Lx-T a correlation between the break 
time T a and the X-ray luminosity Lx measured at T a calibrated (using SNIa) from a sample of Swift GRBs. To reduce the 
uncertainties on cosmographic parameters, since methods like the MCMC are based on an algorithm that moves randomly in 
the parameter space, we a priori imposed some constraints on the series expansions of H 2 (z) and c!l{z), requiring that the 
most general and obvious constraint is the positivity requirement: 

• d L (z) > , 

• H 2 (z) > , 

applied on all our redshift ranges. We perform the same Monte Carlo Markov Chain calculations to evaluate the likelihoods, 
firstly considering the SNIa HD, the BAO and H(z) data sets, or the GRBs- Amati HD, the BAO and H(z) datasets separately 
(Cosmography I), and then constructing an overall data set joining them together (Cosmography II). Instead, we have decided 
not to include the Lx - T a Gamma Ray Bursts Hubble diagram in the overall analysis, performed in Cosmography II, by virtue 
of the peculiar features of this recently discovered correlation, which should be further investigated in order to get stronger 
confidence and confirmation. Either one of the cosmographic analysis are implemented using z (z-Cosmography I /II /III) and 
y = (y-Cosmography I /II /III) series, and, in all the cases, the results are largely compatible. Our MCMC method allowed 
us to obtain constraints on parameter estimation, in particular for higher order cosmographic parameters (the jerk and the 
snap). It turns out that the deceleration parameter clearly confirms the present acceleration phase; both the estimation of the 
jerk and the DE parameters, reflect the possibility of a deviation from the ACDM cosmological model. In particular from the 
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Figure 17. Upper Panel Dependence of w on the j/-redshift in the equation of state described by the parametrization (see Eg. 184 | l for 
different values of the Q m , as described in the label. Bottom Panel Dependence of w on the y-redshift in the same equation of state 
for different values of the cosmographic parameters <jo and jo. The filled region corresponds to the allowed behaviour for the equation of 
state, when the cosmographic parameters are varying within the 2<r region of confidence. f2 m is fixed and set to the value Cl m = 0.22. 




<7o c/o 

Figure 18. Contour plots for the parameter wi(go, jo)i an d for different and fixed values of f2 m , in the case of the CPL parametrization. 
The plots correspond to the values Q m = 0.4, Q m = 0.32 (left and right side-upper panel), Q m = 0.22 and Q m = 0.12 (left and right 
side-bottom panel). The red lines define the 2a region of confidence for the parameters goiJOi as constrained by Cosmography II. 
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Figure 19. Contour plot for the parameter too(<70i^m) in the case of the CPL parametrization. The red lines define the 2cr region of 
confidence for the parameters qo,Q m , as constrained by Cosmography II. 
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Figure 20. Relative residuals between the CPL exactdistance modulus and fifth order series expansions of the same quantity, when 
varying the cosmographic parameters qo, Jo, so, and Zo for the CPL parametrization considered above. The black thick line corresponds 
to the best fit values of the cosmographic values, while the red and green lines define the 3<r region of confidence. The grey triangles 
indicate the observational relative errors for the SNela and GRBs data. 

Cosmography II (which combines the SNIa HD, the Amati Gamma Ray Bursts HD with BAO and H(z) data sets we obtain 
for the parameter jo £ (0.104, 0.92) at 2 a of confidence, and jo € (0.1, 1.4) at 3 a of confidence. We finally investigate the 
implications of our results for dark energy: in particular here we focus on the parametrization of the dark energy equation of 
state (EOS), and we compare the cosmographic and the EOS series. Our analysis indicates that the dark energy equation of 
state is evolving for all the parametrizations we considered; moreover the q(z) reconstruction, allowed by our cosmographic 
analysis, permits a transient acceleration. In a forthcoming paper we are going to compare such indications with a direct 
and full reconstruction of the EOS from the observational data. We showed that the current data sets are not yet able to 
discriminate among these alternative scenarios: the selection of a really high redshift standard rulers is what would really 
improve the knowledge of the expansion history of our universe. 
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